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We present a detailed analysis of strongly driven spontaneous four-wave mixing in a lossy in¬ 
tegrated microring resonator side-coupled to a channel waveguide. A nonperturbative, analytic 
solution within the undepleted pump approximation is developed for a cw pump input of arbitrary 
intensity. In the strongly driven regime self- and cross-phase modulation, as well as multi-pair gen¬ 
eration, lead to a rich variety of power-dependent effects; the results are markedly different than in 
the low power limit. The photon pair generation rate, single photon spectrum, and joint spectral 
intensity (JSI) distribution are calculated. Splitting of the generated single photon spectrum into 
a doublet structure associated with both pump detuning and cross-phase modulation is predicted, 
as well as substantial narrowing of the generated signal and idler bandwidths associated with the 
onset of optical parametric oscillation at intermediate powers. Both the correlated and uncorrelated 
contributions to the JSI are calculated, and for sufficient powers the uncorrelated part of the JSI 
is found to form a quadruplet structure. The pump detuning is found to play a crucial role in all 
of these phenomena, and a critical detuning is identified which divides the system behaviour into 
distinct regimes, as well as an optimal detuning strategy which preserves many of the low-power 
characteristics of the generated photons for arbitrary input power. 

PACS numbers: 42.65.Lm, 42.65.Wi, 42.50.Dv 


I. INTRODUCTION 

Integrated optical microresonators continue to develop 
as a promising platform for generating, controlling, and 
measuring quantum states of light [iHini- Advances in 
fabricating such chip-based structures is enabling the 
construction of micron-scale optical ring resonators with 
reported quality factors Q of over one million [1 [nj. 
By exploiting the nonlinear optical response of the ring 
medium, combined with the massive enhancement of in¬ 
traring pump intensity made possible by the large Q val¬ 
ues of the resonator, a wide variety of nonlinear optical 
phenomena can be realized using relatively modest input 
powers. Entangled photon pair generation in silicon mi¬ 
crorings has been demonstrated using mere /iW of pump 
power 1311 ], and optical parametric oscillation in a sil¬ 
icon nitride microring has been observed using a pump 
power of only 50 mW [12]. Arrays of coupled silicon mi¬ 
crorings have also been investigated as a potential source 
of heralded single photons m- 

Such high-Q microrings are ideal for investigations of 
strongly driven nonlinear optical effects. Depending on 
the application, these effects can be undesirable or highly 
sought after: multi-pair production from an entangled 
photon pair source contaminates the sought after en¬ 
ergy correlation, whereas optical parametric oscillation 
(OPO) only arises in the strongly driven regime. Theo¬ 
retical studies of quantum nonlinear optics in integrated 
microresonators have typically treated the nonlinearity 
perturbatively [HHa], which limits calculations to quan¬ 
tities relating to a single generated photon pair. 

Recently we presented a general theoretical treatment 
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of photon pair generation arising from spontaneous four- 
wave mixing (SFWM) in microring resonators, fully ac¬ 
counting for the quantum effects of scattering losses 
within the resonator El- As our focus was on the ef¬ 
fects of such losses, we limited ourselves to the low power 
regime in which a perturbative solution of the relevant 
equations of motion provides an adequate description of 
the pair generation process. In this work we extend our 
analysis to a more strongly driven regime, where pertur¬ 
bative strategies are inadequate and competing nonlin¬ 
ear effects, including self-phase modulation (SPM) and 
cross-phase modulation (XPM), become important. We 
restrict ourselves to one of the most common pump states 
used in experiment, that of a coherent, narrowband con¬ 
tinuous wave pump beam, for which a nonperturbative, 
analytic solution to the semiclassical equations of motion 
is achievable within the undepleted pump approximation. 
This approximation limits us to pump intensities below 
the onset of OPO, the threshold for which is clearly indi¬ 
cated by our equations; OPO in such structures will be 
the subject of a later communication. Even below the 
OPO threshold, the subtle interplay between the various 
nonlinear terms that couple the ring modes, as well as 
the effects of multiple photon pair generation, give rise 
to a rich variety of nonlinear optical phenomena that are 
accessible by varying only two input parameters, namely 
the pump intensity and detuning. 

In Sec. [TT|we begin by assembling the relevant Hamil¬ 
tonian and field operators for the ring-channel system. 
In Sec. |III| a brief review and summary of our earlier El 
theoretical framework is presented, wherein the system’s 
dynamics are reduced to a set of coupled ordinary dif¬ 
ferential equations for the ring operators alone. Steady 
state solutions for the pump mode, incorporating the ef¬ 
fects of SPM and scattering losses, are developed in Sec. 
|IV[ and the stability of those solutions is studied. The 
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equations of motion for the signal and idler modes are 
then solved in Sec. [V| enabling the calculation of physi¬ 
cal quantities including the photon generation rate, single 
photon power spectrum, and joint spectral intensity dis¬ 
tribution. For each of these measurable quantities the 
corresponding predictions at low and high pump powers 
are compared, and we identify a set of experimental fea¬ 
tures, or “smoking guns,” that distinguish the qualitative 
behaviour at high pump powers from that at low pump 
powers. 


II. HAMILTONIAN AND FIELDS 


We consider an integrated microring resonator side- 
coupled to a channel waveguide, as illustrated in Fig. 
We assume the ring size and quality factor Q have been 
chosen such that the ring accommodates individual res¬ 
onant modes which are well separated in frequency; that 
is, we are in the high finesse limit. While a simple gener¬ 
alization of our framework can be used to treat arbitrarily 
many ring resonances, we restrict our model for the time 
being to contain only three ring modes of interest. The 
full system Hamiltonian can then be written m as 

H — -^channel “1“ -faring -f^coupling -f^bath: (1) 

wherein i^channei refers to the channel fields, i^ring to the 
ring modes, i^coupUng to the coupling between the chan¬ 
nel and ring, and i^bath to any modes into which ring 
photons may be lost, as well the coupling of those modes 
to the ring modes. Introducing channel fields ' 0 j(^), the 
channel Hamiltonian is 


-f^channel — ^ ^ 

J L 
ihvj 


huj j dz ipl{z)'ipj{z) 

d^jiz) 


dz 


dz 


-Ipjiz) - H.c. 


, ( 2 ) 


where the fields satisfy the usual commutation relations 


['ipj{z),ipj>{z')] = 0 , 


^j{z)Ar{z') 


S(z - z')Sjj>. 


( 3 ) 


The index J G {P, S', /} runs over three fields of interest, 
respectively labelled P, S and I for pump, signal and 
idler, with corresponding reference frequencies ujj and 
propagation speeds vj. Each field contains frequency 
components centred at cjj, taken to be the resonant fre¬ 
quency of the corresponding ring mode, and ranges over 
a bandwidth that does not overlap with those of other 
fields but involving excitation over sufficiently long 
distances that the Dirac 6 function in (© is a good ap¬ 
proximation m- By allowing the fields to have different 
propagation speeds we include the possibility of group ve¬ 
locity dispersion between the different channel fields. The 
Hamiltonian © does assume group velocity dispersion 


within each channel field is negligible, but it is straight¬ 
forward to include arbitrary dispersion. The spatial co¬ 
ordinate 2: ranges from ^ — 00 to ^ = +00 with the 

coupling to the ring assumed to take place at a single 
point z = 0 . Within this point coupling approximation 
the coupling Hamiltonian becomes 

-f^coupling ~ E (^7J&jV’j(0) + H.c.) , (4) 

J 

in which we have introduced ring-channel coupling co¬ 
efficients 7 j, as well as discrete ring mode annihilation 
operators bj. In addition to the physical channel, to 
simulate scattering losses in the ring we include an ex¬ 
tra “phantom channel” into which ring photons can be 
lost. The phantom channel similarly accommodates three 
fields (j)j{z) with respective propagation speeds uj and 
coupling coefficients jij^ and is represented as H'bath by a 
channel and coupling Hamiltonian identical to those for 
the physical channel la. 

The Hamiltonian for the ring modes can be written as 

-faring — E hujjb^jbj -I- i^NL, (5) 

J 

where includes all the nonlinearity in the system. 
Since the fields will be most intense within the ring res¬ 
onator, we neglect channel nonlinearities and take ii^NL to 
contain only ring mode operators. In this work we con¬ 
sider effects arising from the third-order nonlinear sus¬ 
ceptibility in the ring, taking 

^NL = (^A.bpbpb\b\ + H.c.^ + Tirjb'^pb'^pbpbp 

+ + b\b'^pbibp^ . ( 6 ) 

The first term is responsible for SFWM, in which two 
pump photons are converted to a signal and idler pho¬ 
ton pair. The second leads to SPM of the pump, while 
the latter two are responsible for XPM between the pump 
and signal and idler modes. It is safe to neglect SPM and 
XPM terms that involve only the signal and idler modes, 
since the power in those modes will be small compared 
to that in the pump mode. While we focus in this work 
on SFWM involving a single pump mode, it is straight¬ 
forward to incorporate multiple pump modes into our 
model. The nonlinear coupling coefficients A, r] and ( 
are not independent, as they arise from the same nonlin¬ 
ear susceptibility, but we formally leave them arbitrary 
for the time being so that the effects of each term in i^NL 
can more easily be identified. Obtaining expressions for 
these constants depends on the approximations used to 
derive the nonlinear sector of the ring Hamiltonian. We 
present our derivation of this Hamiltonian and the asso¬ 
ciated constants A, r] and ( in Appendix \A\ arriving at 
an estimate of 

_ h(^%cn 2 . 

IL V ring 

with T] = Kj 2 and C = 2 A. In this expression 77.2 refers 
to the nonlinear refractive index of the ring material, n 
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^ fj. Phantom Channel 



FIG. 1. Integrated ring-channel system geometry with la¬ 
belled ring modes and incoming and outgoing outgoing chan¬ 
nel fields. Photons generated in the ring may exit to the 
physical channel or be lost to the upper effective “phantom 
channel”. 


to the linear refractive index, and luring to the volume 
of the ring mode. For the silicon nitride rings used in 
typical experiments m, with 77.2 ~ 2.4 X 10-i9m2/w 
[ 22 ] this yields A ^ 10 Hz. For typical silicon rings [5], 
with 77.2 ~ 2.7 X 10“^^m^/W [23] this calculation predicts 
A - 10^ Hz. 


III. EQUATIONS OF MOTION 

The Heisenberg equations of motion for the field op¬ 
erators and (j)j{z^t) and the ring operators hj(t) 

follow from the Hamiltonian 0 , and can be simplified 
by the introduction of auxiliary quantities [T4]: here we 
summarize the results. 

The equations of motion for the channel fields are 

= -i'fjbj{t)6{z), ( 8 ) 

with similar expressions obeyed by the phantom channel 
fields <pj{z^ t). Note that the solutions to these equations 
contain a discontinuity at 2 : = 0 , which is a consequence 
of our point-coupling assumption. To avoid explicitly 
dealing with this discontinuity, it is helpful to introduce 
formal channel fields which we identify as those fields 


which are incoming and outgoing with respect to the 
coupling point. We define the incoming field 
by 

i’j<{z,t) ='tpj{z,t) for z < 0 , (9) 

and extend it to ^ > 0 by requiring everywhere that it 
satisfy the homogeneous version of 0 , 

V’J<(2:,i) = 0. (10) 

This confers a false future on ' 0 j<( 2 ),t), corresponding 
to the free evolution of the incoming field without any 
coupling to the ring. We similarly define the outgoing 
field ' 0 j>( 2 :,t) by taking 

= i’j{z,t) forz> 0 , ( 11 ) 

and demanding that for all ^ 

f d d \ 

( ^+*wjj tAj>(2:,t) =0, (12) 

giving 2 pj^{z,t) a false past to the left of the coupling 
point. By an identical procedure we may define the in¬ 
coming and outgoing phantom channel fields 0 j<( 2 :,t) 
and (/)j>( 2 ), t). Since we will primarily be concerned with 
the properties of the photons generated in the ring, which 
exit to one of the channels and propagate to positive z, 
all calculations involving the ring’s output will be carried 
out on the outgoing fields '0j> (z, t). Our goal is therefore 
to construct an explicit solution for these fields in terms 
of the incoming fields ' 0 j<( 2 ), t). Indeed, since these fields 
freely propagate, the field at large positive 2 ) (where any 
measurements on the generated photons would occur) is 
entirely determined by the outgoing field at 2 : = 0 , 

iPj{z,t) = - ^/vj) forz>0. (13) 

It therefore suffices to construct a solution for '0j>(O,t), 
which can be very simply related to the incoming field 
'0j<(O,t) and the corresponding ring operator bj{t) [T4| 
via 

V’J>(0,i) = V’J<(0,i) - (14) 

vj 

For each operator Oj{t) it will be convenient to define 
the corresponding slowly-varying barred operator Oj{t), 

Oj{t) = (15) 

In terms of these quantities and the incoming and out¬ 
going fields, the equations for the ring mode annihilation 
operators bj{t) are found to satisfy 
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+ Tp + 2iriPp{t)bp{t)^ bp{t) = -ijp'ipp^{0,t) - i/Up^p<(0,i) - 2zA*6p(t)6s(t)6/(t)e“*^™s*, (16a) 

+rs +iC6p(l)6p(i)^ 6s(i) = -i7sV’s<(0)^) - il^*s^s<i^^t) “ iA6p(i)6p(i)6|(i)e*^"“®*, (16b) 

+r7 + iC6Ui)6p(i)^ 6/(i) = -i'yi'tpi^{0,t) - iii*j^j^{0,t) - iAbp{t)bp{t)bl{t)e"^"^^\ (16c) 


where we have introduced the ring mode detuning 

^ring ~ 2CiJp, (1^) 

as well as the total effective linewidths Fj, 

Tj = Fj + Mj, (18) 


where Fj and Mj denote the damping rates associated 
with the physical channel and phantom channel cou¬ 
plings, respectively: 


rj = 


Mj = 


\ij? 

2vj ’ 

M, 

2uj 


(19) 


These total damping rates can be simply related to the 
quality factors Qj of the resonator modes; for example, 
for the pump resonance 


Qp = 


Fp’ 


( 20 ) 


which yields Qp r\j 10^ for a ring with Fp = 1 GHz 
given a pump with wavelength A = 1550 nm. The cou¬ 
pled set of driven, damped ordinary differential equations 
( p^ fully describes the nonlinear dynamics of the ring- 
channel system. Combined with the channel transforma¬ 
tion a solution to this system of equations permits 
the calculation of any measurable quantities on the out¬ 
going photons in the channel. 

It is important to note at this stage that our treat¬ 
ment neglects the effect of ring heating due to the large 
circulating pump power present in the ring. Such ther¬ 
mal effects are routinely observed in experimental inves¬ 
tigations of microring systems, and typically manifest as 
an effective power-dependent drift in the resonant fre¬ 
quencies of the ring as it undergoes thermal expansion 
[12]. For slowly varying and cw pumps a simple way 
to account for this is through the addition of a pump 
photon number-dependent correction to each resonance. 
Our model already incorporates a similar effect: SPM 
and XPM of each mode are represented by precisely 
such terms. The inclusion of thermal resonance drift can 
therefore be modelled by altering the coefficients r] and 
C in Eqs. (16), which would be replaced by effective con¬ 
stants ?7eff and 


'^eff — '^thermal 5 

Ceff ~ C Cthermal- 


While T] and C are both positive, r/thermai and Cthermal 
would be negative, since as the ring expands the res¬ 
onant frequencies are typically lowered [24]. Depending 
on the relative magnitude of the thermal drift coefficients 
compared to the SPM and XPM strengths, in some cir¬ 
cumstances ?7eff and Ceff laay become negative. While for 
the remainder of this work we neglect thermal drift of 
the ring resonances, so that /^thermal = Cthermal = 0, we 
emphasize that our conclusions do not depend sensitively 
on this assumption unless otherwise stated. 


IV. STEADY STATE PUMP SOLUTION 

The set of coupled equations ( [T^ treats both the pump 
and signal and idler modes quantum mechanically, retain¬ 
ing the operator nature of bj{t) for each J. While this is 
necessary if one wishes to fully account for the nonclassi- 
cal properties of the pump mode, in typical experiments 
um the system is pumped by a coherent laser beam or 
pulse. In such situations the initial pump state is de¬ 
scribed by setting each incoming pump mode to a coher¬ 
ent state. The pump field can then be well approximated 
by its expectation value, which is a classical function of 
time. To implement this semiclassical approximation we 
take 


bp{t) I3p{t) = {bp{t)). 


( 22 ) 


In addition to treating the pump classically, we also im¬ 
plement the undepleted pump approximation. In the 
equation for the ring pump mode (16a) the term involv¬ 
ing Opbsbj accounts for the effect on the pump mode 
when a signal-idler photon pair is produced. Neglect¬ 
ing such effects, we drop this term and instead take the 
semiclassical pump amplitude l3p{t) to satisfy 


^ +rp + 2iri\l3p{tp'^l3p{t) 
= -ij*p{'ipp<{0,t)), 


(23) 


in which we have assumed (0p^(O,t)) = 0, so that there 
is no incoming pump energy in the phantom channel. 
Note that while this approximation amounts to neglect¬ 
ing pump depletion due to photon pair generation, linear 
pump losses are still accounted for in our model, as ev¬ 
idenced by the presence of the damping term Fp in Eq. 

d^). 








5 


In this work we consider the case of a continuous wave 
(cw) pump beam injected in to the channel, so that 




IP 


(24) 


where Ap is the detuning of the injected pump from the 
ring pump resonance, and p is a constant related to the 
input pump power Pin in the channel at the coupling 
point via 


p = 


^FpPir: 

Tiujp 


(25) 


In steady state, after the ring pump mode has come to 
equilibrium with the channels, we expect there to be a 
constant average number of pump photons Np in the 
ring, where 


Np = lim \Pp{tp. 


t^OO 


(26) 


Defining /3p(t) = e*^^^/3p(t), from Eq. (23) we have 


— + Fp + i{2r]\Pp{tp - Ap) ) j3p{t) = -ip. (27) 


It is not difficult to show that Np will be constant only 
when I3p{t) has both constant amplitude and constant 
phase, so that dl3p{t)/dt = 0. Setting this time derivative 
to zero in the above equation and taking the modulus 
squared of the result, we find that in steady state Np 
must be a root of the cubic equation 


Cp{Np) = 0, 


(28) 


where 


Cp{Np) = 

- 47?ApA| + (Fp + A%)Np - \pf. (29) 

In the absence of SPM (when 77 0 , or when the input 

power is very small), Np is related to the incoming power 
by a simple linear function. 


Np = ^ 

rl 




Al 


(30) 


The presence of SPM, however, complicates the task of 
determining Np as a function of |pp for a given detuning 
Ap and nonlinearity 77 . The cubic equation ( [^ has in 
general as many as three real, positive roots. Further¬ 
more, only some of these may correspond to stable solu¬ 
tions of (23). Before solving for the roots of Cp{Np), we 
first derive a set of criteria to assess the stability of any 
such solution. 


To determine whether or not a given root of (29) is 


stable, we conduct an analysis similar to that of Hoff, 
Nielsen and Andersen [25]. For a given constant solution 
/3p^ to (27), we define the fluctuation amplitude 6l3p{t) 


/3p(i) = /?(?) +<5/?p(t). 


(31) 


Keeping terms up to first order in d^p^ the equations of 
motion satisfied by 5pp{t) and dpp{t) can be written as 


dt \pi^p{t)) \^/^p(^)) ’ 

where F is the 2x2 matrix given by 


(32) 


F = 

^-Vp - i{ApNp - Ap) 


1 2 




2777 


/j(0)* 


—2iri /3p^ 

—F p + 7(477Ap — Ap) 


(33) 


For a given solution to be stable, we require the real 
part of both eigenvalues of F to be negative, so that 
the fluctuation term SPp{t) will decay with time. These 
eigenvalues are 


f± 


-Fp± 


- {pNp - Ap)2. (34) 


Now, Re(/_) < 0 automatically; demanding that 

Re(/p) < 0 yields the condition 


—2 


Ar]^N^P - {Ar]Np - ApY < F 


p- 


(35) 


Solving this inequality, we And that any solution Np for 
Eq. ( j^ corresponds to a stable solution of Eq. (27) 
if jApJl^s below a “critical detuning”, |Ap| < AcriticaU 
where 


^critical 


= Vsri 


(36) 


When |Ap| > Acriticab ^ solution Np of (28) corresponds 
to a stable solution of (27) if and only if Np lies outside 
a certain interval, Np ypV_, Ap), where 


N± = ^ 
^ 377 


± 2 V 


A2 

critical 


(37) 


Having established the stability criteria for a given Ap, 
we return to the task of finding real, positive roots of (28). 


While analytic expressions for the roots exist, it is more 
instructive to use indirect arguments to study their na¬ 
ture. Taking the derivative of Cp with respect to Ap, 
we And that dCpjdNp = 0 at Ap = A±. Thus, for sub- 
critical detunings (|Ap| < Acriticai), where the Ap are 
not purely real, there are no local extrema - it is easy to 
show that a graph of Cp{Np) is monotonically increasing 
and intersects the Ap axis only once, leading to a sin¬ 
gle real, positive root Ap which corresponds to a stable 
solution. On the other hand, for supercritical detunings 
(|Ap| > Acriticai): iFc functiou Cp{Np) goes through a 
local maximum at A_ and minimum at Ap. The number 
of times Cp{Np) intersects the Ap axis is then deter¬ 
mined by the power parameter |pp; varying it translates 
the graph of Cp{Np) vertically. If Cp{N_) > 0 and 
(7p(Ap) < 0, the graph of the function intersects the Ap 
axis three times, indicating the existence of three real, 
positive values of Ap. The outer two correspond to stable 
solutions, while the inner root is unstable. These multi¬ 
ple cases are illustrated in Fig. in which Ap is plotted 
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FIG. 2. (Colour online) Steady state average photon number 
in the ring pump mode as a function of channel input power 
with Fp = 1 GHz, ?7 = 1 Hz for (a) zero detuning, (b) sub- 
critical detuning, (c) supercritical detuning. Red and blue 
curves indicate stable solutions, green unstable. The dashed 
line respresents choice of optimal detuning to maximize Np 
at each input power (Ap = A^^{Np) = 2r]Np). 


as a function of input power for various detunings. For 
values of |Ap| above the critical detuning of \/3 F p there 
exists a region of optical bistability, in which two stable 
equilibrium average pump photon numbers for a given 
input power are permitted, a phenomenon that has been 
observed experimentally in microring systems [24]. In 
this region the two stable solutions are separated by an 
unstable (and therefore physically inaccessible) range of 
Np. Also plotted in this figure is the case of “optimal 
detuning”, Ap = A^p^{Np), in which Ap is not taken to 
be fixed, but chosen to exactly cancel the effect of SPM 
as Pin is increased, 

Ap = A°P‘(Vp) = 27? Ap, (38) 

which restores the simple linear relationship between Np 


and |p|^ 


Ap = ^. (39) 

Tp 

This behaviour is indicated by the dashed line in Fig. 
[^ which corresponds to a stable pump solution for all 
input powers, always lies on or above the fixed-detuning 
curves, and at each input power corresponds to the choice 
of detuning that maximizes Np. 


V. SIGNAL AND IDLER DYNAMICS 

Having developed the steady state pump solution, we 
return to the signal and idler equations of motion. We 
first develop an exact solution to these equations, valid 
for a cw pump of arbitrary intensity, and then use this 
solution to calculate the photon pair generation rate, as 
well as the one- and two-photon spectra of the generated 
photons. 


A. Exact solution 


We begin by writing the equations (16b-16c) for the 


signal and idler ring operators in the presence of a classi¬ 
cally described cw pump that leads to a ring pump am¬ 
plitude of the form ^p{t) = where pp is a 

constant. Letting hx{t) = for x = S, I we 

obtain 

where M is the 2x2 coupling matrix defined by 


M = 

-^s - ^(CI/^pP - ^p) 

r; 


-7A/?p 


iA/3p _r, + i(c|/3p|2-Ap)) 


and D{t) the driving term responsible for quantum fluc¬ 
tuations from the physical and phantom channels, 


( 7e-*^p*(7/7/)k(0,t) +M/^l<(0,t)) j 


In obtaining M we have assumed the ring resonances are 
equally spaced, so that Aring = ujs ojj — 2ujp = 0; 
the pump detuning, however, is left arbitrary. Previ¬ 
ously El we employed a perturbative approach in the 
frequency domain to solve these equations, while ignor¬ 
ing the effects of SPM and XPM. While this provides an 
adequate description of the pair generation process for 
low pump powers, a nonperturbative strategy is needed 
to treat the strongly driven case. In the cw regime, where 
M is time-independent, this coupled set of linear ordinary 
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differential equations can be solved exactly in the time 
domain for arbitrary pump intensities by taking 

(pS) " / (43) 

—oo 

where the 2x2 matrix Green function G{t, t') is given 
by 


G{t, t') 


g/;, MM” ^ gM.(t-t') 

(gD{t,t') gAit,t')\ 


(44) 


For simplicity we henceforth assume the ring-channel 
coupling constants and propagation speeds for each mode 
are the same, vj = uj = v and /ij = /i so Fj = F 

for each J. The matrix elements go and qa are then given 
by 

X ^cosh[p(t — t')] — sinh[p(t — 

and 

_2 

gA{t,t’) = sinhipit - t% (46) 

in which we have introduced the dynamical parameter p, 
p = - {CNp - Ap)2. (47) 


Depending on the pump photon number Np and detun¬ 
ing Ap, p may be either purely real, purely imaginary, 
or exactly zero. Indeed, as will become clear in the fol¬ 
lowing sections, p serves as an important parameter in 
characterizing the system’s behaviour. 

With explicit solutions written down for the ring op¬ 
erators 6j(t), we can make use of the incoming-outgoing 
channel field relation (14) to determine '0j>(O,t). We 
find for the signal 


and 

qsi{t,t')= (50) 

2 * 

- i')e-»'^p(*+*')e-rO-*') smh[p{t - t')]. 

vp 

The phantom channel response functions are related to 
these via 

Pss{t,t') = ^{qss{t,t') (51) 

and 

Psi{t,t') = -qsi{t,t')- (52) 

7 

Similar response functions pix{t,t') and qix{t^t') can be 
introduced for the idler fields, which, due to our assump¬ 
tion of equal coupling coefficients and propagation speeds 
for the signal and idler fields, are identical to those for the 
signal: pis = Psi, Pii = Pss, Qis = Qsi and qn = qsS‘ 


B. Photon generation rate 


Armed with explicit expressions for the outgoing fields 
'ipgy and we can calculate any measurable quantity 
related to the generated signal and idler photon pairs. 
Of particular interest is the photon pair generation rate, 
one of the primary figures of merit used in assessing the 
practical utility of the ring-channel system. The steady 
state outgoing flux of signal photons Js into the physical 
channel can be calculated via 

Js = lim 

t^OO 

= lim ^ f P\qsi{t,t'p. (53) 

t^OO i J 

Computing the integral, we find 


Js 


2TIS?Nl 


^2 


r -p 




(54) 


t/)5>(0,t) = 


(48) 



qss{t, t')ips< (0, t') + Pss{t, t')(j>s< (0, t') 


+ qsi{t,t')^p\Ao,t') +psi{t,t')(j)\AQA) , 


where we have introduced the temporal response func¬ 
tions qxx' for fho physical channel and Pxx' for 
the phantom channel: 

qssit,t') = 5{t-1') 

- 'Aqu _ p\^-(T+iAp)(t-t') ( 43 ) 

V 

X [cosh[p(t — t')] — - — sinh[p(t — t')], 

P 


The nature of the scaling of Js with pump photon num¬ 
ber Np depends intimately on the character of p, the 
behaviour of which as a function of Np for various de¬ 
tunings is illustrated in Fig. Recalling ( [T7| ), we find 
that p is real when Np G [Ap/^A, Ap/A], with p = 0 at 
the endpoints of this interval, and imaginary otherwise. 
For low enough Ap, when p ^ i\ Ap\, Js scales quadrat- 
ically with the number of pump photons Np. Since in 
this regime Np is directly proportional to the channel 
input power Pin, the overall scaling of Js with Pin re¬ 
mains quadratic, in agreement with experiment [5]. As 
the pump power increases, however, the scaling of Js is 
affected by several separate power-dependent processes. 

First, for a fixed detuning Ap, the SPM-induced drift 
of the pump resonance slows the scaling of Np with chan¬ 
nel input power Pin, as demonstrated in Fig. Second, 












FIG. 3. (Colour online) Real and imaginary parts of p as a 
function of iVp for zero, subcritical and supercritical detun¬ 
ings, indicating that p is always either purely real or purely 
imaginary. The red line indicates p = F. The transition 
between p being purely imaginary and purely real occurs at 
the points Np = Ap/3A and Np = Ap/A. For supercrit¬ 
ical detunings, there exist points where p = F (represented 
on the plot by red diamonds), indicating the onset of OPO 
behaviour. The nonlinear parameter A taken as A = 10 Hz. 



Channel input power [nnW] 

FIG. 4. (Colour online) Photon pair generation rate as a 
function of channel input power for various detunings. The 
dashed curve indicates the optimal detuning case. System 
parameters for this plot are p = 1 Hz, F = 1 GHz. 


XPM between the pump, signal and idler modes effec¬ 
tively shifts the resonance lines of the signal and idler 
modes, compromising the resonance enhancement of the 
pair generation process. Finally, for supercritical detun¬ 
ings |Ap| > Acriticai, P ^ r when Np ^ N±, where 
N± are the same two photon numbers that define the 
stability of the pump solution (37). In that limit the 
photon flux Js formally diverges. This unphysical pre¬ 
diction corresponds to the onset of optical parametric 


oscillation[TTl[T 2 ]. As this threshold is approached, stim¬ 
ulated emission leads to photon pairs being generated 
faster than the rate at which they are removed from the 
ring, preventing the system from reaching a steady state 
within our model. We are prevented from treating this 
case by our assumption of an undepleted pump. Our 
results are expected to be valid when the intraring con¬ 
version efficiency Faring is much less than unity; this ef¬ 
ficiency, defined as the ratio between the steady state 
signal (or idler) and pump photon numbers, can be ex¬ 
pressed as 


F - 

^ring - 


A^ATp 


—2 


r -p 




(55) 


While in future work we intend to investigate the OPO 
regime and the associated effects of pump depletion, for 
the time being we restrict ourselves to regimes where 
Faring ^ 1 ; in all examples presented below this inequal¬ 
ity is satisfied. 

Perhaps most remarkable is the special regime of opti¬ 
mal detuning, wherein Ap is chosen to maximize Np at 
each channel input power, Ap = A'^^{Np) as defined in 
Eq. (38). For this choice p = 0 identically for all A’p, 


p = - iCNp - Ap f 


= JiPNl - (2ANp - 2-Np 




= 0 . 


(56) 


The photon pair flux then maintains its quadratic scaling 
with both Np and channel input power Pin over its entire 
domain: 


Js 


8f3a2 

{fwjpff^ 


P- 


(57) 


This cancellation between the effects arising from pho¬ 
ton pair generation, XPM, and the SPM-dependent de¬ 
tuning strategy A'^^{Np) arises from the simple fixed 
relationship between the associated nonlinear coupling 
strengths A, 77 and Crucial for this phenomenon is 
that the strength of the photon pair generation process 
scale quadratically with the pump photon number Np. 
This cancellation effect would therefore not be possible 
using, for example, spontaneous parametric downcon- 
version, the strength of which would scale linearly with 
Np. The presence of thermal resonance drift would not 
compromise the existence of an A/'p-dependent detuning 
strategy that yields p = 0 over all Np, though such a 
strategy would no longer correspond to that which also 
linearizes and maximizes the relationship between Np 
and channel input power. 

The photon pair generation rate as a function of chan¬ 
nel input power is plotted in Fig. |^for various values of 
Ap alongside this optimal detuning case. For lower pow¬ 
ers, when SPM and XPM are negligible, a pump beam 
with Ap = 0 gives the best scaling of Jp. For inter¬ 
mediate powers the detuning may be tweaked to combat 
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SPM and XPM in order to maximize while for high 
powers the optimal detuning strategy of Ap = A^/^{Np) 
beats any fixed subcritical detuning. The behaviour of 
these curves suggests a simple experiment to identify the 
presence of nonperturbative, strongly driven effects: one 
could simply measure the outgoing signal or idler power 
as a function of pump input power for a set of fixed, 
subcritical pump detunings. For fixed nonzero detun¬ 
ings Ap < Acriticab strougly driven effects are indicated 
by the presence of a global maximum of generated sig¬ 
nal power at intermediate pump input power, followed by 
decreasing signal power approaching an asymptotic value 
of 

2r 

hujs lim Js(Np) = hujs^‘ (58) 

Np^oo 3 

For critically coupled ring systems F ^ r/2 [14], so the 
asymptotic signal power can be related to the total ef¬ 
fective ring linewidth F as simply If thermal 

detuning of the ring resonances is included this asymp¬ 
totic power will be different; however, the qualitative be¬ 
haviour of the signal power as a function of pump power 
will be unchanged. 



FIG. 5. (Colour online) Origin of splitting in the signal and 
idler lineshapes. Green, blue, and red curves indicate pump, 
signal and idler resonances, respectively, and could respresent 
the enhancement factor m that would characterize the ratio 
of the intensity in the ring to the incident channel intensity 
in a linear experiment. The dashed green line indicates a 
pump detuned by Ap, which leads to a generated pair having 
either its signal or idler photon detuned by ~ 2Ap from the 
corresponding resonance, as indicated by the signal and idler 
pairs connected by dashed black lines. The presence of pairs 
from both cases leads to a doublet structure for both the 
signal and idler lineshapes. 


C. Single photon spectrum 

Another physical quantity of interest is the spectral 
lineshape of signal and idler photons that are emitted 
from the ring. For low power cw pumps these single 
photon spectra typically exhibit a Lorentzian lineshape 
UMM with a characteristic width determined by the 
total effective linewidths of the microring cavity reso¬ 
nances. As we now demonstrate, these spectral charac¬ 
teristics are significantly different in the strongly pumped 
regime. We develop results for the signal field spectrum; 
the idler field will have identical properties. 

We define the power spectrum m for the signal chan¬ 
nel field as 


vs{(^s',t)= (59) 

t+T/2 

t-T/2 

where the first-order temporal coherence function 
is defined by 






(60) 


In writing (59) we have introduced the relative frequency 
co-ordinate which corresponds to a frequency offset 
from the ring reference uJs- The physical frequency uJs 
associated with uj. is therefore 


CJs = UJs 


(61) 


In the remaining sections we adopt this notation of low¬ 
ercase subscripts for frequency offsets: uJs for the signal, 
ujp for the pump and uji for the idler. 

Evaluating ( [bOj ) and setting ti = t, ^2 = t-hTwe obtain 

g''^'>{t,t + T) = ( 62 ) 

pcosh[p|T|] +rsinh[p|T|] 

'p - y 

which is independent of t, depending only on the relative 
time difference r, as would be expected for a cw pump. 
Taking the Fourier transform, we arrive at an expression 
for the lineshape, 

vs(ps) = ( 63 ) 

AVVh^Nl 

v^|r -p + i{LOs - Ap)|2|r + P + i{aj, - Ap)|2 ’ 

with an identical equation for the idler lineshape 
This expression takes the form of a product of two 
Lorentzians. We consider first subcritical detunings. 
When p is imaginary, these Lorentzians have identi¬ 
cal characteristic widths = F and are centred on 
UJs = Ap ± |p|. For low powers, when p i|Ap| the 
spectrum is therefore peaked at cjs = 0 and ujg = 2Ap, 
in agreement with the perturbative calculation [T4| . This 
splitting is easily understood as a consequence of the 
tradeoff between energy conservation and resonance en¬ 
hancement of the pair generation process. As illustrated 
in Fig. when a photon pair is produced with a de¬ 
tuned pump, either the signal photon or idler photon in 
a pair, but not both, can be generated within a ring reso¬ 
nance; energy conservation then requires the other to be 
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FIG. 6. (Colour online) Spectral lineshape i^si^s) scaled to unit maximum vs. ring pump photon number Np. For each plot 
we take F = 1 GHz and A = 10 Hz. (a) Ap = 0, (b) Ap = O.SAcriticai ~ 1.4 GHz, (c) Ap = 1.5Acriticai ~ 2.6 GHz, and (d) 
Ap = Ap^^(Ap) = 2r}Np. The origin of the frequency axis corresponds to the ring resonance at us- 


generated with a frequency that lies away from its corre¬ 
sponding resonance. This is seen in the Np 0 limit of 
Fig. [^b). At sufficient pump photon number Np a sim¬ 
ilar splitting can arise from the effective XPM-induced 
detuning of the signal and idler ring resonances even for 
a pump with Ap = 0, as seen in Fig. [^a) for large Np. 
When Ap = 0 the lineshape begins as a singly-peaked 
Lorentzian, eventually splitting to a doublet structure 
when p becomes imaginary as a consequence of XPM. 

For nonzero Ap, as Np increases, XPM effectively 
counters the pump detuning and the extent of this split¬ 
ting is reduced as \p\ decreases, eventually vanishing 
when Np = Ap/3A. If Np is increased further, p be¬ 
comes real and ceases to contribute to spectral split¬ 
ting, resulting instead in an effective correction to the 
linewidth. The lineshape then takes the form of a prod¬ 
uct of two Lorentzians, both centred on ujg = Ap, with 
respective widths 5uj± = F ± p. As p becomes compara¬ 
ble to F the smaller of these two widths becomes domi¬ 
nant, leading to a lineshape with overall effective width 


~ F — p. For subcritical detunings, as demonstrated 
at large Np in Fig. [^b), the spectral splitting is then 
resumed as p once again becomes imaginary. 

For supercritical detunings, as the threshold for opti¬ 
cal parametric oscillation is approached p ^ F and the 
bandwidth of the emitted signal and idler photons be¬ 
comes arbitrarily narrow, as seen in Fig. gc). This fol- 
lows from our idealization of the pump as an indefinitely 
coherent cw beam; in actual experiments the bandwidth 
of the generated photons will become comparable to that 
of the pump, a phenomenon that has been observed in 
strongly pumped experiments on silicon nitride micror¬ 
ings [12]. 

Finally, as shown in Fig. id), in the special case of 
optimal detuning when Ap = A^^(A/'p), so that p = 0, 
the lineshape remains peaked at a single A’p-dependent 
frequency for each Np^ with unchanging characteristic 
width = F, precisely mimicking the low-power result 
at zero detuning. 

Experimentally, measuring the signal or idler lineshape 
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as a function of input power for a nonzero, subcritical 
detuning as in Fig. [^b) would reveal the richness of the 
strongly driven regime, and illustrate the behaviour of 
the p paramater, which incorporates the effects of both 
XPM and pair generation. 


D. Joint spectral intensity 

To assess the degree of spectral correlation between the 
signal and idler modes, it is instructive to study the joint 
spectral intensity distribution of the generated photon 
pairs. While it is straightforward to define this quantity 
for a system driven by a train of weak pump pulses, in 
which multi-pair generation can be neglected, it is a more 
subtle task to craft a sensible measure of spectral correla¬ 
tion in the strongly driven cw regime. In particular, there 
is no single function that characterizes a joint probabil¬ 
ity amplitude of signal and idler photons, since in general 
there will be far more than two photons in the quantum 
state of the signal and idler modes. Furthermore, even 
for weak cw pumps, if one introduces outgoing channel 
annihilation operators cj{ujj) via 

I (64) 

and naively calculates expectation values of the form 
(c^(cjs)c|(ci;i)c 5 (ci;s)c 7 (ci;i)), the idealization of a zero- 
bandwidth cw pump leads to ill-defined expressions in¬ 
volving the square of Dirac 6 distributions. 

To resolve these difficulties, in Appendixwe develop 
a model of a typical experiment used to characterize the 
JSI for weakly driven systems, in which the coincidence 
count rate of signal and idler photons at respective fre¬ 
quencies ujs and uji is measured. We the extend the def¬ 
inition of the JSI to strongly driven systems by defining 
the JSI to equal the calculated outcome of such an ex¬ 
periment for arbitrary input power. This definition re¬ 
duces to the usual result for single-pair output states, 
and serves as a sensible measure of spectral correlation 
between the signal and idler fields. This coincidence rate 
can be written as 





X (V’L(^l)d>(^2)V’S>(i3)^/’/>(i4)) , 

where T(t) is the Fourier transform of a transmission 
function T{uj) that resolves the frequencies of the signal 
and idler photons prior to detection, 

m = I ( 66 ) 

and St is the temporal resolution of the coincidence 
counter. In this expression the spatial dependence of the 


field operators 2 pjy{z,t) has been suppressed; the signal 
and idler arms of the JSI measurement are assumed to 
occur at balanced distances from the ring-channel cou¬ 
pling point. 

The four-time expectation value 

is found to natu- 

rally split into two parts, 

= (67) 


where 


A{ti,t2) = J dt' 


( 68 ) 


The function is precisely the first-order coherence 
function defined in Eq. (60) used to calculate the single 
photon spectrum, 

(69) 

P 


p cosh[/9|f3 - |] + r sinh[p|f3 - ti |] 




r -p 

The A{ti,t 2 ) term, after computing the integrals, is given 
by 


. ^2 


A{ti,t2) = 


(70) 


[ai sinh[p|t2 - |] + a2 cosh[p|t2 - h \ 

where the constants ai and a 2 are defined by 

_ XNp-Ap- 

ai = p — % -^-1, 


a2 = F — i^Np — Ap). 


(71) 


The JSI can therefore be expressed as the sum of corre¬ 
lated and uncorrelated terms, 

— d-corrip^s •) T -^uncorr (ws,a;i), (72) 

where 

dcovvip^S'i^i) — (d^) 

St f f 

J dui J dV2A{vi,V2)f{Us-Vl)f{uJi-V2) 

and 


r ip^s 5 ^i) — 


St 

-iy2)f{0Js - Vl)f{uJs - V 2 ) 


(74) 


/*'■/' 

J du[ j du'^ -i' 2 )n^i - - ^ 2 ) 
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As indicated by their labels, /uncorr can be expressed as 
a separable product of functions of oos and uji , while Icorr 
cannot. Each is expressed as a convolution of the Fourier 
transforms A{ui^U 2 ) and of the A(tiA 2 ) and 

functions, 





dt2 




^iy2t2 



and similarly for ^ 2 )^ with the transmission filter 

function T{ui)T{^ 2 ). The A and functions are deter¬ 
mined by the dynamics of the signal and idler modes in 
the ring, while their convolution with the T functions re¬ 
flects the frequency averaging that arises from the finite 
resolution of a realistic JSI measurement. 

Computing the Fourier transform, we find for A 


A{ui,U2) = V^K^Np5{ui + z/2 - 2Ap) 


1 -i 


CNp-Af 

p 


1 + i 


CNp-Ar 

P 


(iAz^ — r + p){—iAiy — F + p) {iAu — F — p){—iAu — V — p) 


(76) 


with Av = (z/i — z^ 2)/2. The term in parentheses multiplying the S function varies on the scale of F. Assuming that 
the measurement frequency resolution (^cjtrans is much narrower than this, the slowly varying term can be pulled out 
of the integrals in ( [7^ , leaving 

/corr(Ws,Wi) « 5tT^A'^Np[D{uJa - Ap,UJi - Ap)]^ (77) 

1 1 I i C^p-^p 2 

X ^ -|- ' P 

(i{uji - Ap) - r + ~p){-i(u}i - Ap) - r + p) (i(uji - Ap) - T- p){-i{uji - Ap) - T -p) 


where 


D{us,Ui) = 

^ I I + ^ 2 ) f{ 


CJ 


S 


Ui)f{Ui 



The function D{ujs^00i) can be interpreted as the 
“smoothed” version of the Dirac 6{uJs +co’^) distribution, 
and arises from the finite bandwidth of the JSI measure¬ 
ment scheme; D{ujs — Ap^uji — Ap) is sharply peaked 
and uniform along the energy-conserving antidiagonal 
line Co’s + cji — 2Ap = 0 with characteristic width Jco’trans 
(the measurement resolution) in the direction orthogonal 
to that line. 

Finally, taking the Fourier transform of we 

find 


g^^\i'i,-U2) = 6{vi -V 2 ) 


(79) 


X 


JEFA^A^I, 


IF - p + i(z.i - Ap) p |r + p + - Ap)p 


As with A, apart from the 5 function this is slowly vary¬ 
ing compared to the measurement resolution; the term 
multiplying the 5 function can be pulled out of the inte¬ 
gral in Eq. (74). The uncorrelated contribution to the 
JSI /uncorr IS therefore well approximated by 


r 5 ^i) 


6t 

2tt 




T^sYs)T^lYi), (80) 


where i^sYs) and i 2 i{uJi) are precisely the single-photon 
lineshape functions given by Eq. (63) as derived in the 


previous section. The uncorrelated part of the JSI is 
therefore proportional to the simple product the signal 
and idler lineshapes. 


For low power cw pumps, wherein multi-pair gener¬ 
ation is insignificant, the uncorrelated part of the JSI 
/uncorr IS negligible and /corr dominates. The JSI then 
takes the form of a narrow antidiagonal line correspond¬ 
ing to the energy-conserving condition cJs+cj^ —2Ap = 0. 
For Ap = 0, the line is singly peaked, as illustrated in 
Fig. I^a). For nonzero Ap the line is distributed among 
two peaks separated by 2Ap, as evident in Fig. IJb), 
consistent with the single photon spectrum derived in the 
previous section. At higher powers, such as in Fig. [^c), 
this splitting can also arise from XPM-induced signal and 
idler detuning even for a pump with Ap = 0. When the 
splitting is due to XPM-induced signal and idler detun¬ 
ing, the JSI remains centred on the unperturbed ring 
resonances. On the other hand, when pump detuning is 
responsible for the splitting, the JSI is translated by Ap 
along both frequency axes. 

In Fig. I^the uncorrelated contribution /uncorr to the 
JSI is plotted for the same pump parameters as in Fig. 

The weight of the uncorrelated contribution is ex¬ 
tremely small compared to the correlated contribution 
at low powers, as indicated by the scales in Figs, [^and 

but grows to an appreciable level at high powers. For 
Ap = 0, as in Fig. |^a), at low Np the uncorrelated 
part of the JSI displays a single peak centred at the ori¬ 
gin. In the regimes that give rise to split lineshapes, as 
illustrated in Fig. [^b) and|^c), the uncorrelated contri¬ 
bution takes the form of four distinct peaks, symmetri¬ 
cally placed about the centre of the overall distribution. 
Two of these peaks lie on the antidiagonal, overlapping 
with the correlated contribution. The remaining two lie 
on the diagonal, and would therefore appear to violate 
energy conservation if assumed to correspond to signal 
and idler photons that originated from the same pair. 
It is therefore natural to interpret these peaks as corre- 
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spending to signal and idler photons that are detected 
from separate pairs. As these uncorrelated, “non-energy 
conserving” peaks are well separated from the correlated 
part of the JSI, they are uncontaminated by the corre¬ 
lated contribution to the JSI. The properties of photon 
pairs detected in these peaks would therefore be expected 
to differ from those detected in the antidiagonal peaks. 
We intend to investigate such properties in future work. 

The form of the JSI depends qualitatively on whether 
p is imaginary or real, a behaviour we saw earlier in the 
single photon spectrum. When p is imaginary, and thus 
contributes to the frequency terms in the denominators of 
Eqs. 0 and (), a splitting in the JSI appears. When p is 
real, and acts as an effective correction to the linewidth T, 
the JSI is localized to a single line arising from the over¬ 
lap of Jeorr with a single peak in /uncorr- For sufficiently 
detuned pumps, in the regime of real p the uncorrelated 
contribution can be large enough to be visible on the JSI 
plot without exaggeration or scaling. Indeed, for super¬ 
critical detunings |Ap| > AcriticaU as the OPO threshold 
is approached and p ^ T the uncorrelated contribution 
vastly dominates over the correlated contribution, as seen 
in Fig. This is an expected consequence of the rapid 
growth in the photon pair generation rate in this regime 
- multiple photon pairs are generated in sufficiently large 
quantities that joint detection of a signal and idler photon 
originating from the same energy-conserving pair is un¬ 
likely relative to the probability of detecting a signal and 
idler photon which originated from separate pairs and 
thus obey no relationship in energy. Another effect seen 
as p ^ T is the narrowing of the entire JSI distribution to 
a small point-like peak centred on {ug^oji) = (Ap, Ap). 
Within our idealization of a zero-bandwidth cw pump the 
area of this point would be limited only by the frequency 
resolution of the JSI measurement scheme, though in ac¬ 
tual experiments the finite pump bandwidth would serve 
as a fundamental lower bound for the overall extent of 
the JSI. 

Perhaps the most definitive experimental indication of 
strongly driven effects lies in the top-right and bottom- 
left uncorrelated peaks of the JSI distribution for a de¬ 
tuned pump as in Fig. |^b). For sufficient detunings 
these peaks are well separated from the antidiagonal and 
thus easily distinguished from the correlated part of the 
JSI. For low powers, wherein only one photon pair is 
generated in the ring at any given time, they would be 
entirely absent from the measured JSI. As the power in¬ 
creases, any non-spur ions coincidence detection of pho¬ 
tons in these regions indicates multi-pair generation, as 
photons generated in those peaks do not conserve energy 
and therefore must be associated with separate, indepen¬ 
dently produced pairs. 


VI. CONCLUSION 

We have investigated the strongly driven regime of 
spontaneous four-wave mixing in microring resonators for 
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FIG. 7. (Colour online) Correlated part Icorr of joint spectral 
intensity distribution, scaled to unit maximum, of signal and 
idler photon pairs for a pump with (a) Ap = 0,Ap = 10, 
(b) Ap = 3r, Ap = 10, and |c) Ap = 0, Ap = 2 x 10^ 
Ring parameters are taken as F = 1 GHz and A = 10 Hz. 
The splitting evident in (b) arises from the pump detuning, 
whereas in (c) the XPM-induced detuning of the signal and 
idler ring modes is responsible. 
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FIG. 8. (Colour online) Uncorrelated part /uncorr of joint 
spectral intensity distribution, scaled to unit maximum, of 
signal and idler photon pairs. Pump parameters are (a) Ap = 
0, Ap = 10, (b) Ap = 3r, Ap = 10, and (c) Ap = 0,iVp = 
2 X 10®. Ring parameters are taken as F = 1 GHz and A = 10 
Hz. 
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FIG. 9. (Colour online) Joint spectral intensity distribution, 
scaled to unit maximum, of signal and idler photon pairs for 
(a) Ap = 1.5Acriticai with Np = 9.8 x 10^ (90% of OPO 
threshold) and (b) Ap = 1.5Acriticai with Ap = 1 x 10® (95% 
of OPO threshold). Ring parameters are taken as F = 1 GHz 
and A = 10 Hz. 


a cw pump input. A nonperturbative, exact analytic so¬ 
lution to the semiclassical equations of motion within the 
undepleted pump approximation was developed, which 
permits the calculation of any physical quantity related 
to the outgoing signal and idler fields while fully taking 
into account intraring scattering losses. The effects of 
self- and cross- phase modulation, as well as multi-pair 
generation, were found to drastically alter the nature of 
the photon pair generation process at high powers. A 
critical pump detuning of Acriticai = T, where F is the 
total effective linewidth of the ring resonances, was found 
to divide the behaviour of the system into two regimes. 
For supercritically detuned pumps, a region of optical 
bistability of the pump mode is predicted, and a thresh¬ 
old emerges for optical parametric oscillation of the signal 
and idler modes. Pump power-dependent splitting of the 
generated signal and idler photon spectra was uncovered. 
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arising from both pump detuning and cross-phase modu¬ 
lation. In certain intermediate-power regimes, dramatic 
narrowing of the spectral linewidth of generated signal 
and idler photons associated with the approach to opti¬ 
cal parametric oscillation was found. The joint spectral 
intensity distribution (JSI) was analysed, and found to 
consist of separate uncorrelated and correlated contribu¬ 
tions. The correlated contribution is negligible at low 
powers, but becomes significant as multi-pair generation 
becomes appreciable at higher powers. In the regime of 
spectral splitting, the uncorrelated part of the JSI dis¬ 
plays an intriguing quadruplet of peaks, two of which are 
well separated from the correlated part. An optimal de¬ 
tuning strategy was derived in which the pump detuning 
is chosen to exactly cancel the effect of self-phase mod¬ 
ulation at each input power, maximizing the intraring 
pump intensity. By detuning the pump in this manner 
the effects of both spectral splitting and bandwidth re¬ 
duction are eliminated, and the photon pair generation 
rate continues to scale quadratically with the pump input 
even for arbitrarily high powers. 

Three simple experimental tests of our predictions in 
the strongly driven regime were proposed: 

1. For fixed subcritical nonzero detunings the photon 
pair generation rate as a function of input pump 
power is predicted to have a local maximum at 
intermediate powers, followed by a decreasing ap¬ 
proach to an asymptotic level at high powers. 

2. The single photon spectra of the outgoing signal 
and idler fields are predicted to show spectral split¬ 
ting proportional to the pump detuning at low 
powers, followed by a regime of a singly peaked 
spectrum with pump power-dependent narrowing 
of bandwidth at intermediate powers, and finally 
resuming a doublet structure at high powers. 

3. The presence of two non-energy conserving peaks 
lying on the diagonal of the joint spectral inten¬ 
sity distribution, which are a consequence of multi¬ 
pair generation, is predicted to occur for sufficiently 
large pump powers. 

Our analysis was restricted to cw pump inputs; study¬ 
ing how these strongly driven phenomena are altered for 
short pulses requires a numerical approach. Addition¬ 
ally, a slightly more sophisticated solution is required to 
fully study the regime of optical parametric oscillation, in 
which the undepleted pump approximation breaks down. 
We intend to extend our techniques to treat these regimes 
future publications. 


Appendix A: Calculation of nonlinear coupling 
constants 

To estimate the nonlinear coupling constants A, r] and 
C, we present in this section a derivation of the nonlinear 
sector JJnl of the ring Hamiltonian. For the moment 


we imagine the ring has been decoupled from both the 
physical and phantom channels, so that it is idealized as 
a perfect, isolated cavity. We expand the electric field 
E(r) and electric displacement field D(r) in the ring in 
terms of discrete ring modes EQ,(r) and DQ,(r) as 

a 

D(r) = y] + H.c., (Al) 

a 

where are the mode frequencies and the associ¬ 
ated annihilation operators. The contribution to the ring 
Hamiltonian arising from the third-order nonlinear sus¬ 
ceptibility can be written [28] as 

JJnl = (A2) 

-T y drr^^\r)D\r)DHr)D'^{r)D\r) 

with implied summation over repeated lowercase Ro¬ 
man indices, where eo is the permittivity of vacuum and 
represents the nonlinear response coefficients. 
Within the rotating wave approximation, only keeping 
relevant terms for the pump, signal and idler modes, we 
obtain 


^NL = 

I 

4eo 

I 

4eo 

I 

4eo 

I 

4eo 

I 

4eo 


(A3) 
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^^^^^Qppisb'pb'pbibs 


Qppppb'^pb'^pbpbp 


4! 


1!1!1!1 !) V 
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where the constants Qijkl are given by 




Qijkl = (A4) 

I dr (qg*(r)(DHr))*(Dyr))*Z)^(r)Dl(r)) . 

As is typically done for dispersive media we take [29] 



egn2(r; LUi)n‘^{r; uj2)n‘^{r; cj3)n2(r; CJ 4 ) ’ 


where n(r; uj) is the linear refractive index of the ring 
medium at frequency cj, and is the frequency- 

dependent nonlinear susceptibility. To evaluate the coef¬ 
ficients Qijkl we introduce co-ordinates for the ring 
and 1(f), such that the volume element 


dr = pdpd(j)dz 


(A6) 
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can be written as 


dr = dl^ = dr±dl^, 


(A7) 


where R is the nominal ring radius and = R(j), which 
varies from 0 to 27tR = L, the nominal ring circumfer¬ 
ence. The co-ordinate is understood as shorthand for 
the pair (p, z). Writing the mode fields EQ,(r) as 


Ec.(r) = 


eairj)e 


ikalfh 


Vl 


(A8) 


where = 27rnalL for integer Ua^ we can simplify i^NL 
to 


^NL = 


(A9) 


3 huip / huis , t 

J2~^y ~2 ^Qsippbsbibpbp 


L2 2 


3 huip huishui , , t ;,t 

-Qppisbpb'pbsbs 


2 2 
2 


-^(^) Q'ppppbl^bybpbp 
- ^ Q'spspbsbpbsbp 

in which the reduced constants Q'jj^l given by 
1 


Q'l 


IJKL 


(AlO) 


VZiZjZkZl 

J dr_Ld4eox|^3)*(rx, 4 )(e}(r_L))*(e^j(rj_))* 

. (All) 

The modes are normalized [2^ such that 

= (A12) 

1 J dr±dl^eon^{r±-,u}a)el{rj_) • ea(r_L)7gp(r_L; w^), 


1 , 


(A13) 


where 7gp(rx; is the ratio of the group and phase ve¬ 
locities of the ring medium at each frequency an d spa tial 
point. However, we display the explicitly in (AlO) so 


that the expression can be used regardless of whether or 
not the ea{T±) are normalized such that = I- 

To estimate the constants Q'uxl approximate the 
ratio ^gp ^ 1 everywhere, and assume the ea{T±) to be of 
uniform magnitude within the ring and vanish elsewhere. 
We take this uniform magnitude to be unity and assume 
that for the modes of interest {2kp — kj — ks)L 1. For 
modes with polarization mainly perpendicular to the ring 
plane, the relevant susceptibility will be = X(3)^ 

dependent of position in the ring; this can be immediately 


generalized to treat other mode polarizations. We then 
have 


eon^A, 


(A14) 


where A is the cross-sectional area of the ring. Taking 
ujs ^ oji ^ Up in the prefactors of (A9), we finally obtain 


^NL ~ (jiAbpbpb^gb\ + H.c.^ -1- hr]bl^b\^bpbp 

+ hC (^b^sb\,bsbp + b\b\^bibp^ , (A15) 


where 


A = 


3^^pX(3) 


(A16) 


Aeon^LA 

and T] = A/2, = 2A. In terms of the more ex¬ 

perimentally accessible nonlinear refractive index n 2 = 
3x(3)/4eocn^ this becomes 


A = 


hupcn2 


(A17) 


in?LA 

which is in line with the results of similar derivations 


Appendix B: An operational definition of the joint 
spectral intensity 

When characterizing a source of entangled photon 
pairs, the joint spectral intensity (JSI) distribution is of¬ 
ten introduced in the low power limit, when the state of 
the signal and idler modes is well approximated by 

I'i/’s/) = Pvacivac) (Bl) 

+ I duJs j duif{us,Ui)al{us)a\{ui)\Ydic), 

where |pvacP < 1 is a constant and the a^j{u) refer to 
the creation operators at frequency u for the signal and 
idler modes [30]. The unsymmetrized and unnormalized 
JSI for the such a state is defined as \f{us^Ui)\‘^, and 
is proportional to the probability density per unit time 
of jointly detecting a signal and idler photon pair with 
respective frequencies Us and Ui. For strongly pumped 
sources, when multiple photon pairs are generated in sig¬ 
nificant quantities so that higher order terms involving 
more than two creation operators appear in the state, 
it is less straightforward to define a single function that 
characterizes the energy relationship between simultane¬ 
ously detected signal and idler photons. Instead, one 
can operationally extend the definition of the JSI to 
strongly pumped sources by calculating for arbitrary in¬ 
put power the outcome of experiments designed to mea¬ 
sure \f{us^Ui)\‘^ in the low power limit. In this section 
we develop such a calculation for a typical measurement 
scheme employed to measure the JSI of photon pairs pro¬ 
duced in a microring resonator. 

We consider a standard experimental setup m to 
measure coincidence rates between signal and idler pho¬ 
tons of particular frequencies as illustrated in Fig. 
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FIG. 10. (Colour online) Schematic of experimental setup for 
measuring the joint spectral intensity distribution. Pump, 
signal and idler {P, S',/} outputs from the ring resonator 
are incident on a filter F that removes the pump compo¬ 
nent from the beam. Signal and idler helds with modes cs 
and Cl are then separated by dichroic beamsplitter DBS, and 
independently hltered by monochromators, which are imple¬ 
mented by frequency-dependent beamsplitters BSl and BS2. 
Each monochromator-beamsplitter transmits in a small win¬ 
dow (^cjtrans about uji and cjs, respectively. Broadband pho¬ 
todetectors D1 and D2 measure the detector modes C 5 ',det and 
c/,det, and are connected to coincidence counter CC to regis¬ 
ter joint detection events within a temporal resolution of 6t. 
Vacuum is input to the empty ports of BSl, BS2 and DBS. 


The signal and idler fields are separated, and each field 
is sent through a separate monochromator set to trans¬ 
mit photons in some small range (^cjtrans about a centre 
frequency uJq for the signal and uJi for the idler. Placed 
after each monochromator are broadband photodetectors 
connected to a coincidence counter to identify simultane¬ 
ously detected signal and idler photons. The transmis¬ 
sion frequencies uJq and uJi are independently controllable, 
and correspond to a single point (or, more accurately, sin¬ 
gle bin) on the JSI plot, which is produced by scanning 
through ujs and uJi and measuring the corresponding co¬ 
incidence rate. The transmission width is chosen to be 
much smaller than the linewidth of the measured pho¬ 
tons, (^cjtrans T, SO that the full 2D spectrum can be 
resolved. 

The monochromators can be simply modelled as 
frequency-dependent beamsplitters. Provided both the 
signal and idler arms of the experiment are balanced, 
the spatial dependence of the fields after the ring can 
be suppressed; all fields in this section are understood to 
be evaluated immediately after the ring-channel coupling 
point. We introduce annihilation operators csioog) and 
ci{uJi) for the ring output fields ' 05 >(t) and as in 

Eq. ( [ 6 ^ . We can then apply the appropriate transforma¬ 
tion to obtain the annihilation operators cj^^eti^j) for the 
fields seen by the detectors placed after the monochro¬ 
mators. For the signal, we have 

OS',det(^) — T(cd CJg)cs (co) + F(uJ ^s)^S,vac(^) -> (B2) 

in which cs,vac refer to the modes on the other input port 
of the monochromator-beamsplitter, into which only vac¬ 
uum is present. The transmission and reflection functions 


T(u;) and R(u;) determine which frequencies are transmit¬ 
ted by the the monochromator. For example, a simple 
filter may be modelled by a transmission function with a 
rectangular frequency profile. 


f(co) 



< 

otherwise. 


<^^trans 

2 


(B3) 


The exact choice of T(cj) is not important for our pur¬ 
poses; for simplicity we only assume T(cj) is a real, suffi¬ 
ciently narrow, symmetric function of cj. The reflection 
function will be irrelevant, though it will satisfy the usual 
restrictions to correctly model a beamsplitter. 

In exactly the same manner, modes c/^det(<^) seen by 
the detectors of the idler arm can be introduced. We can 
then write down the fields measured by each detector in 
the usual way, 

(B4) 

In a typical coincidence measurement the signal detector 
is continuously activated, and detection of a signal pho¬ 
ton at time t is used to trigger the activation of the idler 
detector (which is placed at a small delay relative to the 
signal detector) for a very short time St, so that the idler 
detector samples the idler field during the time interval 
[t — St/2,t -h St/2]. The average rate I{uJs,oJi;t) at time 
t of coincident detection events at oos and uJi of the sig¬ 
nal and idler detectors is given by the standard Glauber 
formula involving the fields at each detector [32] . 


I{uJs,uJi;t) = 

t+T/2 t'+(5t/2 

t J dt' J dt" 


(B5) 


lim 

T^oo 


t-T/2 t'-5tl2 


(V^kdet (^OV^kdet (^0'^7,det i^")) 


(B 6 ) 


In steady state the expectation value depends only on 
time difference \t" — t'\, which in the integrand is at most 
St. Provided the coincidence resolution time St is much 
smaller than the timescale on which the expectation value 
varies (in our case F ), this expression for /(cd^, t) is 
then well approximated by 

(B7) 

V^St {'Ips,det (0'^/,det (0'^S',det (0'^/,det (^) ) • 

Proceeding to expand the detector fields in terms of 
their constituent modes, we find 


I{uJs,uJi;t) = 

v‘^St 


(B 8 ) 


(27r) 




dz/4 




X T*(z/1 - LOs)f *{l>2 - C0i)f{v3 - UJs)f{Vi - Wj) 
X {c\{vi)c\{v2)cs{v3)ci{v4)) . 
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By writing the operators cj{ui) in terms of their respec¬ 
tive parent fields and then carrying out the integration 
over each we arrive at Eq. ( [6^ , where T(t) is the 


Fourier transform of the transmission function T(co’); see 
Eq. 

In obtaining (65) we have again used the fact that in 
steady state the expectation value is invariant with re¬ 
spect to time translations by t in each argument. The 
expression (65) is manifestly real, independent of time, 


and in the limit of small (^cjtrans indeed reproduces the 
single-pair JSI |/(cc;g,when calculated with an ini¬ 
tial state of the form ( |B1| ). 
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